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The effective Lagrangian of a test particle, interacting within an ideal gas, is cal¬ 
culated within the closed time path formalism in the one-loop approximation and 
in the leading order of the particle trajectory. The expansion in the time deriva¬ 
tive, available for slow enough motion, uncovers diffusive forces and decoherence 
in the particle coordinate basis. The master equation, generated by the effective 
Lagrangian, is derived and its consistency is verified for a finite-temperature gas. 


I. INTRODUCTION 

The time-reversal invariance of effective forces is a nontrivial issne: Even if the dynamics 
of a closed system is time-reversal invariant, the effective interaction between a snbsystem 
and the rest, its environment, always breaks the time-reversal invariance dne to the envi¬ 
ronment bonndary conditions in time. A fnrther layer of complications is fonnd in inhnite 
systems where the effective forces may be dissipative. Qnantnm systems raise an additional 
question about the decoherence, another irreversible process, generated by the effective 
interactions. These issues had already attracted much attention, and the emergence of ir¬ 
reversibility has been demonstrated in a physically appealing way, by coarse graining [1]. 
In a similar manner, decoherence can be generated by the large, highly degenerate environ¬ 
ment [2, 3]. Such a general picture naturally leaves open the details about the actual loss 
of information, realized by the use of a restricted set of observables, in a given microscopic 
model. 

The traditional approach to Brownian motion comes from kinetic theory; the quantum 
theory is inferred from a master equation. The master equation, obtained hrst for a particle 
colliding with a gas [2], was shown to describe the decoherence while dissipation has been 
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ignored. Though the resulting divergence for long time can be eliminated without taking 
account of dissipation [4], the physically satisfactory description requires the presence of 
both the dephasing and the dissipation [5, 6]. A systematic, perturbative derivation of the 
master equation, based on the collision cross section of the test particle, taken in the Born 
approximation, has been worked out, too [7-9]. This approach was later improved by going 
beyond the Born approximation [10]. The traditional many-body technique is available, as 
well, to arrive at a master equation [11]. 

Another way to approach this problem is model building. The simplest, exactly solvable 
model consists of a particle, coupled linearly to inhnitely many harmonic oscillators [12-14]. 
The effective dynamics of an open system, the test particle in the present case, can easiest 
be handled within the closed time path (CTP) scheme. This formalism was first introduced 
in quantum held theory [15] and has since then been successfully applied in different areas 
of condensed matter physics [16-18] and particle physics [19]. The harmonic model [21] 
has thoroughly been studied in this scheme [22] and a general form of the master equation 
for the reduced density matrix, containing a memory term, has been derived in the CTP 
formalism [23]. 

Though the Markovian approximation is believed to be applicable at high enough tem¬ 
peratures, the high temperature master equation, derived in Ref. [13], failed to preserve 
the positivity of the density matrix [24-26]. A Markovian master equation of the Lindblad 
form [28], i.e., satisfying the physical requirements for the density matrix, has been found 
for high enough temperatures [29], fitting into the set of minimally invasive generalizations 
[26, 30] of the master equation of Ref. [13]. A physically satisfactory master equation can 
be constructed in an ad hoc manner [31] and by the use of an appropriate ansatz for the T 
matrix in the collisional approach [7], as well. 

The goal of this paper is a systematic derivation of the effective dynamics of a test 
particle, moving in a gas. We shall consider an ideal gas, the only interaction in the model 
being between the test particle and the individual gas particles. The fact that in this model 
dissipative forces occur is not surprising since the leading-order, one-loop contribution to the 
transport coefficients which is equivalent to using Kubo’s formulas [20] can be interpreted as 
the contribution of an ideal gas. The effective Lagrangian of the test particle is derived within 
the framework of the Landau-Ginsburg double expansion, assuming that the interaction with 
the gas generates a small amplitude, slow modihcation of the test particle trajectory. The 
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effective action of a particle which interacts with an ideal gas was calculated some time 
ago by using the traditional effective action approach in imaginary time [32] and by means 
of the CTP formalism [33]. Though our procedure is similar to the one followed in these 
works, the hnal form of the effective action is different. This is because we had to go 
beyond the traditional action formalism to recover dissipative forces and did the separation 
of the conservative and the dissipative forces in an appropriate way. Another element of this 
work is that the decoherence of the coordinate has been monitored and the result has been 
compared with those of the harmonic model [13]. It is found that both the friction forces 
and the decoherence are generated by the same quantity in our one-loop effective theory, 
supporting the view about the common origin of dissipation and decoherence. 

Although the collision-based and the CTP effective action descriptions have different 
starting points they share some common assumptions. In the collisional picture one assumes 
a dilute gas to truncate the hierarchy equations, it requires short range interactions to rely on 
scattering processes, and it needs high enough temperatures to render the effective dynamics 
Markovian. Finally, the Brownian limit is carried out by letting the ratio of the gas particles 
and test particle mass go to zero to arrive at a simpler equation of motion for the reduced 
density matrix. These assumptions are mirrored and partially softened in the effective 
action approach where the effective action is usually constructed within the framework of the 
Landau-Ginsburg double expansion. In fact, the expansion in the amplitude is reminiscent 
of the assumption of weak interactions in a dilute gas. The expansion in the gradient, the 
assumption that the effective interactions are local, corresponds to the use of scattering 
processes in the other scheme. The calculation of the higher order contributions is not 
exceedingly difficult in the effective action, leaving open the way to a systematic inclusion 
of multi-particle correlations and non-Markovian effects. Notice that the light gas particle 
limit is not necessary for the calculation of the effective action. 

The calculation of the order O and O {d'j) effective Lagrangian is presented below. 
The real part displays Newton’s friction force and a mass renormalization whereas the imag¬ 
inary part describes the decoherence of the coordinate. The result simplihes considerably 
when an ideal gas of fermions is considered at vanishing temperature: Both the friction force 
and the strengths of decoherence stem from the same source and are proportional to the 
particle velocity. Furthermore, the decoherence shows a characteristic anisotropy, namely 
it is six time stronger in the direction of the velocity than in the perpendicular plane. The 
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master equation for the density matrix, generated by the effective Lagrangian, is derived 
and its Lindblad form is verihed for hnite temperature fermion gas. It is remarkable that 
one needs the full O {d^) decoherence dynamics to establish this result. 

The paper starts with a short summary of the CTP formalism, presented hrst for a 
closed system in Sec. IIA when it appears as a trivial rewriting of the traditional, transition 
amplitude based formalism of quantum mechanics. The Green’s functions are introduced 
for a single harmonic oscillator in Sec. IIB. The CTP formalism becomes unavoidable for 
open systems when the system-environment interactions are non-conserving and make the 
system dynamics open. The resulting effective theory is briefly introduced in Sec. IIC. Such 
an effective theory is calculated for a harmonic toy model in Section III. The main topics of 
this work, the effective dynamics of a test particle within a gas is embarked in Sec. IV. It 
starts in Section IV A with outlining the perturbative derivation of the effective theory for 
the test particle. The general equations, obtained there, are used in Section IV B to find the 
O (x^), O ((9^) effective Lagrangian, the special case of an ideal gas of fermions at vanishing 
temperatnre being touched upon in Section IV C. The master eqnation, corresponding to 
this Lagrangian, is derived and its Lindblad strnctnre is verihed in Section V. Finally, 
Section VI contains our conclusions. 

II. CTP FORMALISM 

The CTP formalism has been proposed to deal with expectation valnes in the Heisenberg 
representation [15] rather than transition amplitndes of the Schrodinger pictnre. Thongh this 
is already an important extension compared to the traditional, transition amplitnde based 
formalism the real strength of this scheme becomes visible for open systems, interacting 
with their environment. The effective dynamics of such systems, such as dissipation, can 
not be captnred by relying on pure states only and the use of mixed state and the rednced 
density matrix is needed. This is the point where CTP formalism becomes nnavoidable, 
being the only scheme to deal with mixed states in many-body systems and to visnalize the 
elementary processes by borrowing the clarity of the Feynman diagrams. 
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A. Closed system 

Let us suppose that our closed quantum system is described by the coordinate x and a 
time independent Hamiltonian, H. The time evolution operator, 

U{tf, U; j) = T[e-s 4 ' (1) 

where T denotes the time ordering and the source, j(t), a book-keeping device, is written 
in the Heisenberg representation and allows us to recover the expectation value of an x- 
dependent observable, A{x), in the form 

(A(x(t))) = Tr[H(a;(t))p] = Tr[H(a;)H(t, L; 0)piH'f(f, 0)], (2) 

where pi denotes the density matrix of the initial state. The calculation of this quantity is 
facilitated by the use of the generator functional, 

etmil = -J-)], (3) 

where the CTP doublet, j = (j"'", j“), contains two independent sources, introduced for each 
time evolution operator in Eq. (2), and serve to generate the observable and the interaction 
vertices of the perturbation series by the functional derivatives 5/5j^{t), e.g. 

(A(xm = A (4) 

with a = ±. Note that the unitarity of the time evolution assures that either source 
can be used to generate the expectation value for arbitrary tf > t. The path integral 
representation of the generator functional, (3), can be obtained by the slicing of both time 
evolution operators according to the strategy of the traditional path integration for transition 
amplitudes. The result is an integration over the pairs of trajectories, x(t) = {x~^(t), x~(t)), 



where the action, 

S'[a;] = [ dt[L{x~^{t),x~^{t)) — L{x~{t),x~{t))] + S^[x], (6) 

Jti 

contains the Lagrangian of the system, L{x,x), and 

S[x]=i^ f dt[{x+{x)y+ {x-{x)f], 


(7) 
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which keeps the path integral convergent for large x. The boundary conditions in time, 
suppressed in Eq. (5), involve the integration of the initial coordinate with the weight 
p{x^iti), x~iti)) and the hnal condition, x^{tf) = x~(tf), representing the initial state and 
the trace in Eq. (3), respectively. 

To recover a time translation invariant action we have to take the limit ti —oo, 
tf ^ oo. However, this is a rather non-trivial procedure [34] which can be avoided by (i) 
considering a harmonic oscillator, as discussed below, (ii) obtaining its Green’s function 
within the operator formalism, (iii) calculating the inverse propagator, (iv) identifying it 
with the kernel of the harmonic action, and hnally, (v) using the kinetic energy which is 
obtained in such a manner as the kinetic energy for the non-harmonic case. The result is 

= T / dt[{x+{t)f - (a;"(t))2] +-P [ dtdt' ^ (8) 

4 J-oo t-f 

where P denotes the principal part prescription. 

A distinguished feature of the CTP formalism, the redoubling of the degrees of freedom in 
the path integral, follows from the double appearance of the physical state, a bra and a ket, 
in the expectation value, (V’|A|'^). These states move in time in opposite directions, making 
it possible to represent effects which are not time-reversal invariant, and are called chronons 
because both degrees of freedom are needed to describe a general, not necessarily time- 
reversal invariant interaction and thereby to realize the flow of the time as we imagine it [34]. 
The internal time arrows of the two chronons are opposite of each others thus the chronon 
exchange, {x~^,x~) —)■ t{x~^,x~) = {x~,x~^), represents the time-reversal transformation and 
the relation, 

W[j] = -W*[tj], ( 9 ) 

expresses the time-reversal symmetry of the CTP formalism [35]. 

B. A single harmonic oscillator 

The Green’s functions play an important role in solving models and they can easiest be 
introduced for a harmonic oscillator. We consider for this end a harmonic oscillator, dehned 
by the Lagrangian Lq = mox‘^/2 — muj^x^ j^. Its action can be written in the form 

‘5'o[^] = X / dtdt'x(t)D~^(t — t')x(t'), (10) 

^ J — OO 
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and the corresponding generator functional, 


shows that 




^o[j] = ~2 dtdt'j{t)D{t - 


(0|T[a;(t)a;(tO]|0) (0|a;(tOx(t)|0) 

(0|a;(t)a;(t')|0) (0|T[a;(f)a;(t)]|0)’ 


= ihD{t — t'), 


can be interpreted as the CTP Green’s function, with the Feynman propagator in its diagonal 
blocks and the Wightman function in the off-diagonal elements. The Green’s function is 
dehned in the frequency space by the Fourier integral 

D^= [ dte^^^b{t). (13) 


Eq. (12) yields = Gui{oJo)/mo, with 


G.(G) = 




—27ii6{u‘^ — r2^)0(—ca) 


-27ri5{u'^ — r2^)0(a;) 




in a straightforward manner. The inverse of the Green’s function gives, after inserting it 
into Eq. (10), the action (6) in the limit ti —oo, —)■ cx) with (8). 

The initial state in the calculation of the Green’s function (14) is the ground state. If the 
initial state contains excitation then the Green’s function is modihed and the calculation, 
leading to Eq. (14), is to be repeated. For an important family of excited initial mixed 
states, corresponding to thermal equilibrium one hnds 


(^^^(G) —)■ Gij(G) — i27i6{u'^ — r2^)n(a;) 


where n{u) denotes the occupation number. 

The time-reversal symmetry, (9), imposes the relation (u) = (—u))*, which 

implies the characteristic GTP block structure, 

[-f iD\ —D^ + iD{\ 

D=\ , . ( 16 ) 

\ -|- iDl —D^ + iD\ I 

for the Green’s function of any local operator in time in terms of four real functions of the 
time. The expectation value of the coordinate, belonging to a physical source, = ±j. 


{x{t)) = J2^' / dt'D^-'{t-t')3{t'), 



is independent of the choice of a, according to Eq. (4) and the Green’s function of operators, 
obeying linear equation of motion simplihes to 

[ + iD'^ —Df + iD^ 

D=\ ^ , ( 18 ) 

\Df + iD'" —D^ + iD'^ I 

due to the a-independence of the expectation value (17). Due to the symmetry, (u) = 
we have D'^{u) = D"'(—ca), D^{u) = —D^{—uj) and D*(a;) = D*(—ca). The 
combination of the Green’s function blocks, appearing in Eq. (17), gi'''’6s the 

retarded Green’s function whose even and odd components, G"' and G-^, are traditionally 
called the near and the far Green’s functions in electrodynamics. The positive norm of the 
states, contributing to the spectral function, 

iD~^ = ^(0|a;_^|n)(n|a;a;|0), (19) 

n 

and the positive excitation energies make D~~^ = 0 for ca < 0 and D '“(a;) > 0 if ca > 0 and 
induce the relation 

iDl = sign(a;)Dl. (20) 

The generator functional of a harmonic oscillator is actually given by Eq. (11) up to 
a source independent constant, corresponding to the path integral with vanishing external 
source. Due to the unitarity of the time evolution the generating functional lTo[j], defined 
by (3) is vanishing for physically realizable sources, = —j~, lTo[j, —j] = 0. The block 
structure (18) makes jDj = 0 in this case and cancels the source independent constant in 
Wo[j]. 


C. Open system 

The generator functional, given by Eq. (3) can be factorized into the product of the 
transition amplitude and its complex conjugate for a closed system which is in a pure state. 
Thus the GTP scheme has no advantage over the traditional formalism in this case. But this 
changes when an open system is considered. Let us suppose that we observe the coordinate 
X which interacts with its environment, described by the coordinate y, and the dynamics of 
the full, closed system is described by the action S[x,y] = + Se[x,y]. The generator 
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functional (3), 


= / D\x]D\v\e^ 


^S[x+,y+]-^S[x ,y ] + ^Se[x] + ^Se[y\ + ^ J dtj{t)x{t) 


( 21 ) 


assumes the form 


= y i:)[£]e^^<=//[h+s>5.[x]+i/dt£(t)j(t)^ ^'22) 

with Seff[x] = S's[a;’''] — S's[a;“] + Sinfi[x], where the influence functional [36] is given by 


^\Si„fl[x\ ^ / £)J^JgfSe[x+,J/+]-fS'e[x ,y ] + i5. 


(23) 


One hnds a physically better motivated form of the effective action by developing the 
CTP formalism for the density matrix itself, 




(24) 


rather than for its trace, c.f. Eq. (3). This is the open time path scheme which is based on 
the generator functional for the reduced density matrix. 




{n\U{tf,ti;j+)pU\tf,ti;-j )\Xf)^\n), (25) 

n 

where the sum is over an environment basis. The path integral expression. 


p{x^,Xf) = / D[x]D[y]ef^ 


LS[x+,y+]-i:S[x-,y-\ + \SPx\ + iS, 


(26) 


contains integration over chronon trajectories with hnal conditions x{tf) = Xf and y^{tf) = 
y~{tf). The reduced density matrix is given by an effective path integral. 


p(x^,Xf) 




(27) 


and it is important to bear in mind that snch an effective open time path dynamics contains 
the effective action of the CTP formalism. 

To separate the dynamically different terms in S'e// we use the decomposition. 


^eff[x] = 5'i[a;^] — ] + 5*2 [h]. 


(28) 


with S'^S 2 [x]/Sx^Sx~ 7^ 0. Here the action Si describes a closed, conservative dynamics. The 
contribntions to S 2 arise from the presence of several contribution to the sum over the envi¬ 
ronment basis elements in Eq. (25). Thus the couplings between x^ and x~ in S 2 represent 
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the mixed state contributions to the path integral (27), arising from the system-environment 
entanglement and make the effective system dynamics open and nonconservative. 

The imaginary part of the effective action describes two different phenomena. ^82 con¬ 
trols the suppression or the enhancement of the contributions to the density matrix and is 
a function of the separation of the chronon trajectories. In particular, if 9 ^S' 2 [a;+, a;“] 00 

as I —)■ cx) then the density matrix is small for large \x~j — xj\ and the coordinate is 

strongly decohered. The elementary excitations of the effective dynamics are dehned by Si. 
Their life-time is hnite if ‘^Si > 0; however, the unitarity of the time evolution is preserved 
by the mixed terms of the density matrix, owing to S 2 7 ^ 0 . 


III. TOY MODEL 


To prepare the calculation of the effective theory for a test particle in a gas we consider 
hrst a simpler problem within a system of linearly coupled harmonic oscillators [13], defined 
by the action S = Sg + Sg, 

,2 

.2 / / tu 

2 '^ 


^ m .n mUJr. 9 

S, = dt{ —x^ - :Ax^ 


Sg = 




m .2 mw; 2 


2 2 


Un 9nXyn 


(29) 


It is more advantageous to separate the system and environment in a way that the system 
potential reflects clearer the effective dynamics for small x, 


Sg = 


Sg = 


f , 

m 0 1 

/ dt 

2 ■ ( 


mUJn 


E 


9r. 


2mu^ 


X 




m 


-Vn - 


2 / \ 2 

mW„ ( QnX 

Vn + 


moj. 


2 2 

The simplest parametrization of the model is given in terms of the spectral function. 


(30) 


p(^) = 


9 I 


2 mu, 


-(S(w„ - fi). 


(31) 
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A. Effective dynamics 


The usual way of finding the effective system dynamics is based on the Heisenberg equa¬ 
tions of motion, 


j = mx + mulx -t- ^ gnVn, 

n 

0 = myn + mulyn + gnX. 

The elimination of the environment coordinates by their equation of motion yields 


(32) 


yn(l^) = 


9n 


7x{u) 


(33) 


m[{uj + uy — 

where the ie term takes care of the initial conditions. The substitution of this equation into 
the hrst equation of (32) leads to the effective equation of motion, —i{oj) = {D'^{uj))~^x{uj), 
in terms of the retarded Green’s function, 

1 


= 


m[{uj + iey — Uq] — ’ 

containing the retarded self energy. 


(34) 




m[{u + iey — 


(35) 


The CTP formalism reproduces this simple result within the Lagrange formalism. The 
action is of the form 

= ^ J dtdt’[x{t)D~^{t,t')x{t') +y{t)G~^{t,t’)y{t')] - J dtx{t)agy{t), (36) 

where D(t,t') = G(t — t',uo)/m and Gn(t,t') = G(t — t',Un)/m, denote the system and 
environment Green’s functions, respectively and the GTP matrix. 


a = 


1 0 
0 -1 


(37) 


denotes the metric tensor of the simplectic structure of the GTP formalism. A variable in 
a Gaussian integral can be eliminated by its classical equation of motion, so the effective 
action is given by 

1 


Seff[x] = 


dtdt'x{t)[DQ^{t,t') — aTi{t,t')a]x{t'). 


(38) 
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The self energy, S = gGg, can simply be expressed in terms of the spectral function, 

1 

= — / dn2np{n)G^{n), ( 39 ) 

m Jo 

in particular 


K = —Pf = -^7rsign(a;)p(|a;|), = -7rp{\uj\). 

m Jo uj^ — \ ¥ 

The influence functional, given by the self energy piece of the effective action, is 

Sinfi[x\ = J 


(40) 


(41) 


The CTP formalism provides no more details than the traditional solution, discussed at the 
beginning of this section. But in the presence of interactions the CTP scheme is necessary 
in order to generalize Wick’s theorem for open systems and to represent the contributions 
to the perturbation series of the Green’s functions in terms of CTP Feynman graphs. 


B. Effective Lagrangian 

To derive the effective Lagrangian we assume that the self-energy is a real, analytic 
function of iu and write the influence functional after a partial integration as 

1 _ °° (Y f 

Sinfi[x\ = “2 / dtx^{t)dlx^'{t). (42) 

aa ' ^=0 

The CTP block structure, (18), results in the influence Lagrangian, 

Linfi = -]^{xfrx'^ + x'^frx -f x'^t^x - xt^x'^ + xHfJx^), (43) 

where the parametrization = a; ± x^/2 is used and the arrow on 

E"”' = E (44) 

£=0 

is a reminder that this is a differential operator and it acts to the right. The effective theory 
is unitary by construction and this property remains valid if the series (44) is truncated. This 
can be seen in an explicit manner by noting that the Gaussian integral, (22), can be found 
by evaluating the integrand at its saddle point, satisfying the classical equations of motion. 
In the case of a physical external source, j*’" = aj, we have x'^it) = x~(t), or x‘^(t) = 0, 
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which cancels the hnite, O {e^) part of the action (6) and the influence Lagrangian, (43). 
The resulting equation, W[j, —j] = 0, assures that the probability is conserved for arbitrary, 
physical external source. 

In the case of a harmonic system, the expectation value, (x), satisfies the classical equation 
of motion and the variation equation for x'^ at x'^ = 0 is 


m{x) = —{mujQ + E’’)(a;), (45) 

with E*” = S" + E-f. Up to order O {d^) we have a shift of the harmonic system potential. 


Wq —)■ cug — Aca^, with 


= — 


dfl 




n ' 


(46) 


corresponding to the system potential in (30), and a renormalization of the mass, m —)■ 
m + Sm, where 


6 m = 


m 


pin) 


(47) 


Furthermore, a Newtonian friction force, F = —kx, with k = 7rp'(0) is generated. The 
imaginary part of the influence Lagrangian (43) yields a suppression factor, with 

do = '7rp(0) > 0, in the path integral (27) and generates decoherence in the coordinate 
diagonal basis. 


IV. TEST PARTICLE IN A GAS 


Consider now the more realistic system of a point particle of mass m, moving in a potential 
V[x) and interacting with a homogeneous ideal gas in equilibrium via a time-independent 
potential, Uix). The action, S = Sp + Sg + Si, is the sum of the terms 


Sp[x = 


Si[x,^lJ\^lJ] = 


^x‘^{t)-Vixit)) 


dt 

dtd^yilj^{t,y) 
dtd^yU{y - x{t))'il;\t, y)ip{t, y). 


ihdt -\ -A -|- /i 

2 m 




(48) 


A. Effective theory 


The influence functional. 






( 49 ) 
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contains the action of the gas particles, 


Sg[ijj\'ijj] + Si[x,'ijj\'ijj] = / dxdx'%i)\x){F ^{x,x')+ V[x,x']x])%i){x'), 


(50) 


where 


with [37] 




ujq 


uj—Uq-\-ie 

—i27rd{u — Uq) 


0 

1 


^nqi2nS{u — u, 





(51) 


(52) 


and Cq = hq^/2mg, Uq denoting the occupation number. The exchange statics factor is 
^ = +1 for bosons and —1 for fermions. When the Green’s function is calculated in the 

operator representation then the excited state contributions in the initial state are picked 
up when the normal order, aa^ = + 1, is used and they generate the overall factor, 

^Uq, in front of the on-shell propagator of the thermal bath excitations. The temperature is 
supposed to be sufficiently high in the case of a bosonic ideal gas to suppress the condensate. 
The particle-gas interaction is described by the term in the action, where T^^'[x\ = 

The exact influence functional, corresponding to an ideal gas. 


Sinfi[x\ = —ihTrln[F ^ r[a:;]]. 


(53) 


is non-polynomial in the trajectory. By assuming that the test particle-gas interaction is 
sufficiently weak one can use U to organize a perturbation expansion in which the leading, 
O (G^), order effective action is [33] 

Sinfi[x\ = ^jaGaj, (54) 

after canceling the O {U) tadpoles by the introduction of a hctitious, homogeneous, neutral¬ 
izing charge density. The source, = U{y — x'^it)), denotes the perturbation of the 

ideal gas by the test particle, and 

G'^^'^^{xi,X2 ) = i^hnsF'^^'^‘^{xi,X2)F'^^'^^{x2 ,Xi) (55) 

is the density two-point function, given in terms of the Green’s functions = G^g + Gl^^q, 
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^iq — ~ C^-cjq) — G+ G_^q, with [37] 




G 


Loq 




d?k Uk 

(27r)3a; - Uk+q + uJk 
d^k 

j^;^nk{nk+q + OS{(^ 


^k-\-q H“ ^k) •) 


(56) 


where the spin is taken into account by the degeneracy factor, Ug = 2s + 1. Since iG~g = 
—‘2G~q is the spectral function for the excited states of the gas G~^ < 0 and G~g = 0 for 
ca < 0. At hnite temperatures the integrals (56) are analytic functions of the dimensionless 
variables x = uj/\q\vF, y = |q|//cF, G = G{x,y), here vp = fikp/rrig, kp = ^GTT^n/n^, n 
being the density. 

The O {U'^) influence functional, (54), is non-local in time and non-polynomial in the 
trajectory. A simpler structure can be recovered by resorting to the Landau-Ginsburg double 
expansion: We assume that the external potential, V{x), acting on the test particle is 
sufficiently attractive, allowing us to use the displacement, generated by the interactions 
with the gas, |a;|, as a small parameter. The O {x^) approximation of the effective theory 
retains the leading self-energy correction to the test particle and contains dissipative forces 
which can be generated by a linearized equation of motion. Another small parameter is 
generated by restricting our attention to a slowly moving test particle where the kernel 
of the linearized equation of motion can be expanded in the time derivative. Such an 
expansion renders the effective action local in time which is expected to be the case if the 
gas has sufficiently high temperature. When the influence functional (54) is evaluated in 
Fourier space then the dominant contribution comes from the region of u which is around the 
characteristic frequency of the particle motion, u ~ \qx\-, indicating that the particle should 
be much slower than the characteristic velocity of the gas, \x\ <^vp and |a;| = |a;|/|q|ni? 1. 


B. Finite temperature 


We are in the position to work out the effective Lagrangian, corresponding to the O {x^) 
truncation of the influence functional (54), written as 


Sinfl = - 


dfdv I (i+1 it- 

( 2„)4 . 


(57) 
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It yields the influence Lagrangian, 


Linfl 




where the notation x = x{t), x^'^^ = d^x/df^ is applied. We assume that the potential V(x) 
is sufficiently deep to confine the particle close to its initial position, x = 0, and retain the 
leading order, O (x^) part only, 

, f d^q 




a' \'\ 2 /^aa' 
Oq • 


The derivatives in 


A±a; = ^ 


n=l 


2^n! 


(59) 


( 60 ) 


act on the Green’s function, evaluated at vanishing frequency. The block structure (18) 
gives the expression 

roo 


Linfl 


127r2 


poo 

/ dqq^U^ {AAixAix'^Gq^ + AAix{x'^ + A2X^)G, 

Jo 


d\nf 

Oq 


+ [A^x^^ - {x^^ - A 2 xy]tGl^] 


(61) 


after integrating over the direction of the wave vector q and the insertion of Eq. (60) brings 
the influence Lagrangian into the form 


Linfl = ^ ^ g 2 m-l, 2 n- 2 X 


m,n=l 


m,n=l 


{ 2 m—l)^d{ 2 n— 2 ) 


+t Y1 ^ d2m-l,2n-lX^^^-^^ , (62) 

m,n=l m,n=l 


where the effective coupling constants. 


g2m—l,2n—l — 


g2m—l,2n—2 — 


d2m,— l, 2 n—l 


3^222(m+n-3)(2m - l)!(2n - 1)! 
i 

37r222(m+n-4)(2m _ i)\{ 2 n - 2)! 

1 r 

127r222(™+»^-3)(2m - l)!(2n - 1)! Jq 

( _/*oo 


2jj2^{jn+n')—2^jy. 


dqq^U^di 


q 




■ 27 ' 7'2 ;:)2(m+n)—2^_ 
'^0g> 


dqq^U^A 


q ILO 


d2m 2,2n 2 STr222im+n-3) (^2m - 2)\{2n - 2)\ J q 


2tt 2 Pi2(m+n)—4 q— 


0, (63) 


have been introduced. The real terms in the Lagrangian contribute to the equation of 
motion and the even and odd order derivative terms, generated by G” and G-^, represent 
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time-reversal invariant interactions and friction forces, respectively. The imaginary part of 
the effective Lagrangian originates from G* and governs the decoherence. The constants 
(63) introduce the time scales Tm,n = and 

with m -|- n 7 ^ 2 in the effective dynamics. The time-reversal parity of the term in the 
Lagrangian which is multiplied by gm,n is (—1)™+"' therefore Tm,n is an intrinsic time scale 
of the conservative or the dissipative effects for n + m even or odd, respectively. The 
characteristic time scales of the decoherence, are generated by time-reversal invariant 
terms. 

It is important to realize that both the irreversibility and the decoherence stem from 
the same source, namely from the Green’s function component G~. If the time arrow 
of the gas is flipped then G~ changes sign, indicating that G~ represents the dynamical 
breakdown of the time-reversal symmetry, the dynamical origin of dissipative processes [34]. 
Therefore dissipation appears with different time scales in the effective dynamics, some of 
them corresponding to irreversibility, others belonging to decoherence. The common origin 
of irreversibility and decoherence can clearer be seen by noting that irreversibility has a 
genuine quantum manifestation, the leakage of the quantum state into the environment. 
This is characterized by the life-time of the state, given by the imaginary part of the inverse 
Feynman propagator in the frequency space. The self-energy of the test particle is given by 
the contributions a = a' on the right hand side of Eq. (59) and the life-time arises from 
The decoherence is driven by the imaginary part of the terms a = —cr', containing 
The Green’s function G^^^' possesses the block structure (18) hence both the life-time 
and the decoherence are driven by the same Green’s function component, namely by G*; in 
particular they have the same characteristic time scale. 

The higher order derivatives in the equation of motion create difficulties because the 
corresponding initial conditions are unknown [34]. Hence we restrict ourselves to the order 
O ((9^) where the effective Lagrangian, 

= (m„ + - kx^i + \\,k(x‘f + (64) 

contains the parameters 

5m = gi^i, k = —gig, do = do^, d 2 = di,i (65) 

One can prove that the dynamics, induced by such a truncated Lagrangian, is unitary by 
repeating the argument, mentioned above in relation to the Lagrangian (43). 
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The Euler-Lagrange equation for at {x^) = 0, 

m{x) = —(VE(a;)) — k{x) + O (df) + O ((x^)) , (66) 

given by ^Linfi, shows the mass renormalization, m = tuq + 5m, and the presence of a 
non-vanishing Newtonian friction force. Note that the O {k) friction term generates changes 
of the momentum by the time without violating the translation invariance of the influence 
functional. In the calculation of ref. [33] G-^ was used in the place of G” and so no mass 
renormalization is found. The dissipative term of the equation of motion (66) has been 
found also by the collisional scheme, the friction constant being expressed in terms of the T 
matrix [9]. 


C. Electron gas at vanishing temperature 


If the temperature approaches zero in a fermionic ideal gas then the Green’s function, 
Guiq, simplihes and develops singularities. The loop-integrals in eqs. (56) can easily be 
calculated for Ug = Q^kp — |q|) [37] and G~^, the Lindhard fnnction, is fonnd to be ana¬ 
lytical for |a;| +y/2 < 1. Within this domain we have G~ = — kpmgQ{x)x/27111^, yielding 
G^ = —ikpmgx/2nk? and G* = —kpmg\x\/2Titi^. The influence Lagrangian, using the same 
approximation as before, is now of the form 

= —kx'^x + Smx^x + O [df) + O [x^) , (67) 


with the friction constant 

ml ^ 

and 5m, given by (65). Observe that there are no higher order irreversible terms owing to 
D^{x,y) ~ X. 

In calculating ^Linfi the expansion of the exponent on the right-hand side of Eq. (58) 
must be carried out around iu{qx — u), yielding 


Note that Ax, which now starts at the order O can be ignored dne to G* ~ |a;|. The 
block structure, (18), gives 


^Linfl — . 


d?q 

(27r)3 


Ul{qx 


d\2 


G l I 

q{x+^),q q(x - ;^),q 


(70) 
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which can be replaced by 

Ql- — f d q jj2Qi (qx'^'f ('71') 

when the contributions beyond O {x^) are neglected. The integration over the solid angle is 
straightforward and leads to a particular anisotropy, 

^Linfi = i\\x\{xf + Qxf) (72) 

where x^^ was separated into longitudinal and transverse components, xf, xf, respectively, 
dehned by x'^ = x'^ + xf, x^x = 0 and 

/*2/c77' 

^= 64 ^/ 

It is easy to understand the origin of the factor |i| in (72). The dissipative processes 
and the decoherence are generated by an energy exchange with the environment and this is 
always possible if the gas is in thermal equilibrium. For vanishing temperature the gas is 
in its ground state and the particle must possess some kinetic energy to interact with the 
gas. This condition is automatically satished by the friction force, being proportional to the 
velocity. But the decoherence would start in zeroth order if the particle is coupled to the gas 
by its density. The role of the factor |i;| is to suppress the system-environment correlations 
for a particle at rest. The initial state of the toy model is not at the energy minimum which 
explains that the decoherence is present even for a static system. 

The absence of further, higher order derivative dissipative forces and of the decoherence 
for a particle at rest remains valid when the higher powers of the coordinate are retained in 
the effective Lagrangian as long as the gas is eliminated in the one-loop approximation, i.e. 
Eq. (54) applies. 


V. POSITIVITY OF THE DENSITY MATRIX 

When an approximation is applied one has always to check whether the consistency con¬ 
ditions for the solution remain valid. The effective theory for the test particle describes the 
reduced density matrix which has to satisfy a number of conditions, such as the preservation 
of the total probability, the Hermiticity, and the positivity. It has already been mentioned 
at the end of Sec. IIB that the block structure (18) makes the total probability preserved 
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for a harmonic oscillator. It is easy to check that the same block structure of the Green’s 
function is recovered when calculated in the one-loop approximation and it guarantees the 
Hermiticity of the density matrix, too. Thus what is left to check is the positivity. This is 
easier to do by inspecting the equation of motion for the density matrix, generated by the 
effective Lagrangian (64). 


A. Master equation 


The nowhere differentiable nature of the trajectories, dominating the path integrals, 
makes it important to use the mid-point prescription in the case of coordinate and ve¬ 
locity couplings [38]. This applies to the O {k) friction term of the Lagrangian which is 
generated by G^{—u) = —Gk[u) and therefore determined by the symmetric time deriva¬ 
tive, —kxf^{xn+i — Xn-i)/2At for discrete time steps where /„ = f(ti + nAt). The expression 


x'^(xn+i - Xn-i), written as x'^{xn+i - x^) 


^n-l\S^n 

r»d 


rtf»d N I /-tn rf»d 

^n-l) “T 


Xn-ixf^_^, can be 


replaced by xf^Xn+i — XnX'^j^^ + Xnxf^ — Xn-ix'^_^ in the coupling of two subsequent time 
slices in the path integral. This step leads to the regularized Lagrangian, 


LB{Xn+i,Xn) = m 


iS^n+l ^n){S^ 


d 

n+l 


X^ 


rtf* fY*d I rr* rnd rv* rY*d 


-V{x 


n+l J 


- Vix-^,) 


At 




/y*d _ /ytd ^ ^ 


At 


(74) 


We hrst identify the bare Liouville operator, Ub, leading the time evolution over a discrete 
time step At. Its action on the density matrix can be written in the path integral formalism 
as 

IdBPiA) = I + y): ( 75 ) 

the normalization being given in terms of the bare mass, niQ. Working up to O {At) one 
hnds 




2h ® 


X / 2 hitV'^^ I X _l_ yV-|--(yV)^ I p(;r, t). (76) 


1 


It is easy to check that any other distribution of V(x°') and x'^^ between the two consecutive 
time steps in the Lagrangian (74) yields the same results in O (At). The Gaussian integration 
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in Eq. (76) leads to the differential equation 

dtp = +Vp, 

with H = /2m + V{x), the role of the generator T> being played by 


T) d k , id2k 

T^b = - X Va-d H-^ 

2 n m 


<k j_ doh , 3 , mn 

2 2 m? At m 


(77) 

(78) 


and d = do + d 2 k^/m‘^ denoting an effective decoherence strength parameter. 

Let us now assume, for the sake of the argument, that the particle moves in a harmonic 
potential, V{x) = mu‘^x‘^/2, and its action can be recast in the form (10). Since the Green’s 
function, G, appearing in the influence functional (54), possesses the block structure of 
Eq. (18), the generator functional (3) of the particle is vanishing for physically realizable 
sources, = —j~, even if it is interacting with the gas within our approximation. This 
result indicates that the total probability is preserved in our calculation as can be realized 
by the renormalization, 

V = (79) 

At m 

which takes into account the trajectory independent contributions to the effective La- 
grangian, neglected in the derivation. Since these are independent of V{x), the result 
remains valid for arbitrary V{x). The action of the renormalized generator on the density 
matrix can be written in a more useful form. 




The master equation, derived within the T-matrix formalism [5, 7] is of this form except that 
it is lacking the last term. The latter represents an interference between the friction and 
the O {d/) decoherence terms of the effective Lagrangian (64) and is not explicitly present 
in the collision based approach. 

The equations of motion originating from the quadratic Lagrangian (64) are satished by 
the expectation value of the coordinate, m{x) = —muj‘^{x) — k{x). It is easy to hnd the 
equation of motion for the expectation values of the canonical operators for an arbitrary 
potential, V{x), 


m{x) = {p) 

{p) = -(yv{x)) - ^{p), 

III 


( 81 ) 



22 


which lead to the Euler-Lagrange equation of the effective Lagrangian (64), 


m{x) = —(VE(a;)) + k{x). 


B. Positivity 

The density matrix preserves its positivity as long as the the master equation can be 
written in the Lindblad form [27, 28], 

p = p] + i 5^ (Kp, Vj] + pVji) . (83) 

j 

The generator (80) of the time evolution can be written in such a form with the help of the 
operator set {x,p, a, b}, where 

mwox + ip muox + p 

a = _ 6 = _ :, (84) 

\/2mh\uo\ ^/2mh\uo\ 

Uq standing for an arbitrary frequency, 

Vp = ^ [AH, p] + D^{[xp, x] + [x, px\) + Dp{\pp, p] + [p, pp]) 

+Da{[ap, a^] + [a, pa)]) + -Db([hp, 6] + [b, pb]), (85) 


^ ^ d uk f 1 do) 

HD^ = -- + — , 

2 2 V2 mj ' 

~ O-2 d2 — — [ - -1 , 

2m^ [_ uj \2 m / 

Da = sign(a;)^, 

2m 

, N dok , , 

Db = sign(a;) —. (86) 

The resulting equation is a special case of the possible Lindblad-compatible master equations 
which are bilinear in the canonical variables [39]. To assure Dx,Dp > 0, we have to impose 
the inequality 

/ 1 1 \ _ 2md , 

k I 7^ H ) < u) < , , X, (87) 

\2d2 mj k[m + 2d2) 
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which is always possible as long as 


u <2. 




( 88 ) 


m{m + 4 ^ 2 ) 

Note that both decoherence parameters up to 0{d^), do and ^ 2 , are needed to establish 
friction and one needs a hnite temperature for an ideal fermion gas to generate the positive 
density matrix in our approximation. 

We may generate the O {k) surface term in the Lagrangian (74) by a time independent 
basis (gauge) transformation, iIj{x) —)■ Qil){x) = e^^‘^il){x)^ and use the Lagrangian, 

(nf* nr* /-f* ^ /yt^ 

*^nA*^n+l *^n) ^ •^n*^n+l 


^B(*n+l,*n) = m- 


At2 


“^(*n+l) + ^(*n+l) + 2 

The calculation followed above gives the generator 
ik , , ik'^ 


2 At 

do(a;)^+i)^ + ^2 


^ n+i - 

At 


d \ 2 


V'p = 


2mh 


{ppx - xpp) + 


8mh 


x^p — px^) — 


do 


2w?h 


(pV + PP^ - 2ppp) 


~ 2xpx) 


dok , , ?>k 

l^^i^PP - ^PP - ppx + ppx) + ^P, 


with d' = do + dok"^/‘^rn? which can be brought into the Lindblad form (85) by 

AH = -—x^ 

8m 

d' djkf do 

Da = sign(£7)^, 

2m 

Db = sign(A) 


k 

(. do\ 

do — TT— 1 


2u 

\ m J 


2m? 


The coefficients can be chosen positive if 

V <2 


dod- 


00-2 


m{m + 2 ^ 2 ) ’ 
and the corresponding equations of motion are 


m(x) = (p) - t{a;) 

(p) = -{VV^(a^)) 


k \ 


(89) 


(90) 


(91) 


(92) 


( 93 ) 
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If p is eliminated then the equation of motion for x agrees with Eq. (82). Notice, however, 
that the gauge transformation Q leads to the representation p p + kx/2, which triggers 
the breaking of translation invariance by AH. 

VI. CONCLUSIONS 

A simple scheme was presented to calculate the effective Lagrangian of a particle which 
performs a slow, small amplitude motion in a static potential and interacts with an ideal gas. 
The calculation, performed on the one-loop level, is based on the Landau-Ginsburg double 
expansion in the particle trajectory and the time derivative. Irreversibility appears in two 
different manners in the effective theory: The friction arises from the time-reversal odd 
terms of the real part of the effective Lagrangian and the decoherence originates from the 
imaginary part. A Newtonian friction force and a mass renormalization were found in the 
linearized equation of motion, considered in the O {d^) order. Furthermore the decoherence 
of the particle coordinate has been established up to the same order. In the case of a Fermi 
gas at vanishing temperature the higher derivative dissipative forces and the decoherence of 
the particle momentum are absent in the one-loop effective dynamics. 

It is instructive to compare the dissipative forces with the Abraham-Lorentz force, the 
“friction force” of the electromagnetic radiation held. The Newtonian friction force, dis¬ 
cussed in this work in the context of the electron gas can be considered as the radiation 
reaction force, generated by the emission of particle-hole pairs. This friction force allows the 
measurement of the velocity with respect to the environment and it appears if the system 
looses either the Galilean- or the Lorentz-boost invariance due to the coupling to its envi¬ 
ronment or due to the initial conditions. If the electromagnetic held is boost invariant in 
the initial state then the dihusive part of the radiation reaction force must contain at least 
the third time derivative of the charge trajectory. 

In the case of a hnite temperature fermion gas and a sufficiently weak friction the master 
equation, generated by our ehective Lagrangian can be brought in the Lindblad form. It is 
remarkable that both the full O {d^) Lagrangian and a hnite temperature environment are 
needed to establish this result. 

Unfortunately, a detailed comparison of the coefficients occurring in the master equa¬ 
tions obtained by diherent approaches is highly non-trivial due to the diherent assumptions 
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involved. Although the master equation obtained in the Born approximation shows some 
superficial similarity with Eq. (65), important differences exist because the latter is not 
based on the scattering picture. The derivation of the master equation within the harmonic 
toy model is usually based on the spectral representation of the environment oscillators. 
Our master equation (80) follows as soon as the influence Lagrangian (43) is truncated at 
O ((9^), yielding the effective Lagrangian (64). 

The higher order terms of the Landau-Ginsburg expansion can be included in our proce¬ 
dure, while this task represents a real challenge in the collisional approach. A more careful 
treatment of the non-Markovian memory terms arising naturally in our scheme is possible, 
and it should provide a non-trivial, physically motivated example of non-Markovian behavior 
[40]. 
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